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ISO-mERMALELASTOHYDRODYNAMIC LUBRICATION OF POINT CONTACTS 


IV -STARVATION RESULTS* 
by Bernard J. Hamrock and Duncan Dowson^ 

Lewis Research Center 

SUMMARY 

The theory and numerical procedure developed by the authors in an earlier publi- 
cation were used to investfeate the influence of lubricant starvation on minimum film 
thickness. This study of lubricant starvation was performed simply by moving the inlet 
boundary closer to the contact center. From the results it was foimd that, for the range 
of conditions considered, dimensionless inlet distance at the boundary between the fully 
flooded and starved conditions m* can be expressed simply as m* = 1 + 3. 06 (R /b) 
fJ ' > where R is the effective radius of curvature, b is the semiminor axis 

of me contact ellipse, and H „ is the central film thickness for fully flooded condi- 

r .2 1 0. 56 

tions, or m* = 1+3.34 (R^/b) H^. -.J , where is the minimum film 

thickness for fully flooded conditions. That is, for a dimensionless inlet distance m 
less than m*, starvation occurs; and for m ^ m*, a fully flooded condition exists. 
Furthermore, it has been possible to express the central and minimum film thicknesses 
for a starved condition as H„ <, = H „ Rm - l)/(m* - 1)1®' and H . o = H . 

[_(m - l)/(m* - 1)J^‘ . Contour plots of the pressure and the film thickness in and 

around the contact are shown for the fully flooded and starved lubricating conditions. 
From these contour plots it was observed that the pressure spike becomes suppressed 
and the film thickness decreases substantially as starvation increases. 


Presented at Joint Lubrication Conference cosponsored by the American Society of 
Mechanical Engineers and the American Society of Lubrication Engineers, Boston, Mas- 
sachusetts, October 5-7, 1976. 
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INTRODUCTION 


It was not until the late 1960's and early 1970’s that the influence of lubricant star- 
vation upon elastohydrodynamic behavior received serious consideration. Prior to this 
time it was assumed that the inlets were fully flooded . This assumption seemed to be 
entirely reasonable in view of the minute quantities of Ixibricant required to provide an 
adequate film. However, in due course it was recognized that some machine elements 
suffered from lubricant starvation. 

How partial filling of the inlet to an elastohydrodynamic conjunction influences 
pressure and film thickness can readily be explored theoretically by adoptii^ different 
starting points for the inlet pressure boundary. Orcutt and Cheng (ref. 1) appear to 
have been the first to proceed in this way for a specific case corresponding to a partic- 
ular experimental situation. Their results showed that lubricant starvation lessened 
the film thickness. Wolveridge, et al. (ref. 2) used a Grubin approach (ref. 3) in an 
analysis of starved elastohydrodynamic lubricated line contacts. Wedeven, et al. 

(ref. 4) analyzed a starved condition in a ball -on -plate geometry configuration, and 
Castle and Dowson (ref. 5) presented a numerical solution for the starved line -contact 
elastohydrodynamic situation. In both references 2 and 4 the analysis yielded values of 
the proportional reduction in centerline film thickness from the fully flooded condition 
in terms of a dimensionless inlet boimdary parameter. 

Only in recent years have complete solutions of the isothermal elastohydrodynamic 
lubrication (EHL) of point contacts been presented. The analysis requires the simulta- 
neous solution of the elasticity and Reynolds equations. The authors' approach to the 
theoretical solution is presented in two previous publications (refs. 6 and 7). The first 
of these publications (ref. 6) presents an elasticity model in which the conjunction is 
divided into equal rectangvilar areas with a uniform pressure applied over each area. 

The second (ref. 7) gives a complete approach to solving the elastohydrodynamic lubri- 
cation problem for point contacts. 

The most important practical aspect of the EHL point-contact theory (ref. 7) is the 
determination of the minimum film thickness within the contact. That is, the mainte- 
nance of a sufficient fluid film is extremely important to the operation of the machine 
elements in which these thin, continuous, fluid films occur. In a recent report by the 
authors (ref. 8) the fully flooded results obtained from the theory given in references 
6 and 7 are presented. In the results the influence of the ellipticity and dimensionless 
speed, load, and material parameters on minimum film thickness was investigated. 
Thirty -four different cases were used in obtaining the fully flooded minimum-film- 
thickness formula. 

The present report uses the basic theory developed in references 6 and 7 in stu(fyir^ 
the effect of lubricant starvation on pressure and film thickness within the conjunction. 
The objective of this work is to provide a simple expression for the dimensionless inlet 
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boundary distance . This inlet boundary distance defines whether a fully flooded or a 
starved condition exists in the contact. An additional objective is to develop simple ex- 
pressions for the minimum and central film thicknesses under the starvation condition. 
Fifteen cases, in addition to three presented in reference 8, were used in obtaining the 
starvation results. To more fully understand what occtjts in going from a fully flooded 
condition to a lubricant starvation condition, contour plots of pressure and film thickness 
in and around the contact are shown. The theoretical findings are compared with pre- 
viously reported results. 


SYMBOLS 


A 

a 

B 

b 

E 

E’ 


constant defined in eq. (2) 
semimajor axis of contact ellipse 
constant defined in eq. (2) 
semiminor axis of contact ellipse 
modulus of elasticity 


1 - -1 , 1 - ■'I 


E, 


E, 


F normal applied load 

G dimensionless material parameter, E’/Pi„ „c. 

1 At ^ ^3 

Hg P dimensionless film thickness at center of contact for a fully flooded conjunc- 
tion, hjj/Rjj 

H„.„ dimensionless minimum film thickness, h^.„/R^ 

min ’ min' x 

H . -IT dimensionless minimum film thickness for a fully flooded conjimction, 
min, r ^ j 7 

Vin, f/**^x 

H . o dimensionless minimum film thickness for a starved conjunction as obtained 
min,D 

from least-square fit of data, S^x 

h film thickness 

k ellipticity parameter, aA> 

m dimensionless inlet distance 

m* dimensionless inlet distance at boundary between fully flooded and starved con- 

ditions 
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“w 


p 

p 

%v, as 
R 

r 

U 

u 

W 

x, X,xl 

y, Y,Yj 


dimensionless inlet distance boundary as obtained from Wedeven, et al. 
(ref. 4) 

dimensionless pressure, p/E' 
pressure 

asymptotic isoviscous pressure 
effective radius of curvature 
radius of curvature 

dimensionless speed parameter, utiq/E'R^ 
surface velocity in x-direction 
dimensionless load parameter, F/E’R^ 

coordinate systems defined in report 


a pressure -viscosity coefficient 

7?^ atmospheric viscosity 

V Poisson's ratio 

Subscripts: 

A solid A 

B solid B 

x,y coordinate system defined in report 


BOUNDARY BETWEEN FULLY FLOODED AND STARVED CONDITIONS 

Figure 1 shows the computing area in and around the Hertzian contact. In this figure 
the coordinate X is made dimensionless with respect to the semiminor axis b of the 
contact ellipse, and the coordinate Y is made dimensionless with respect to the semi- 
major axis a of the contact ellipse . The ellipticity parameter k is defined as the 
semimajor axis divided by the semiminor axis of the contact ellipse (k = a/b) . Because 
of the way the coordinates X and Y are made dimensionless, the Hertzian contact el- 
lipse becomes a Hertzian circle regardless of the ellipticity parameter. This Hertzian 
contact circle is shown in figure 1 with a radius of unity. The edges of the computing 
area, where the pressure is assumed to be ambient, are also denoted. In this figure the 
dimensionless inlet distance m, which is equal to the dimensionless distance from the 
center of contact to the inlet edge of the computing area, is shown. 
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Lubricant starvation can be studied simply by reducing the dimensionless inlet dis- 
tance. A fully flooded condition is said to exist when the dimensionless inlet distance 
ceases to influence in any significant way the minimum film thickness. When starting 
from a fully flooded condition and decreasing m, the value at which the minimum film 
thickness first starts to change is called the fully flooded - starved boundary and is de- 
noted by m* . Therefore, lubricant starvation was studied by using the basic elastohy- 
drodynamic lubrication point-contact theory developed in reference 7 and observing the 
effect of reducing the dimensionless inlet distance. 

Table I shows how changing the dimensionless inlet distance affects the dimension- 
less minimum film thickness for three groups of dimensionless load and speed param- 
eters. For all the results presented in this report the material parameter G is fixed 
at 4522 and the ellipticity parameter is fixed at 6. In this table it is seen that as the di- 
mensionless inlet distance m decreases the dimensionless minimvim film thickness 
decreases. 

nun 

Table II shows how the three groups of dimensionless speed and load parameters 
affect the location of the dimensionless inlet boundary m* . Also given in this table are 
the corresponding values of the dimensionless central and minimum film thicknesses for 
the fully flooded condition as obtained by interpolation of the numerical values. The di- 
mensionless inlet boundary distance m* shown in table II was obtained by using the data 
from table I when the following equation was satisfied : 


^min, F ' (®min) 


m=m* 


= 0.03 


^min, F 


( 1 ) 


The value of 0.03 is used in equation (1) since it was ascertained that the data in table I 
could only be obtained to an accuracy of 3 percent. 

The general form of the equation that describes how the dimensionless inlet distance 
at the fully flooded - starved boimdary varies is given as 


m* 



( 2 ) 


The right side of equation (2) is similar to the forms of the equation given by Wolveridge , 
et al. (ref. 2) and Wedeven, et al. (ref. 4). By using the data obtained from table I, the 
following equation can be written: 
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2 


-,0.58 


m* = 1 + 3.06 



H 


c,F 


(3) 


A fully flooded condition exists when m ^ m* , and a starved condition exists when 
m < m* . 

If in equation (2) the dimensionless minimum film thickness is used instead of the 
central film thickness, the following equation is obtained: 


m* =1 + 3.34 



-,0.56 


^min,Ij 


(4) 


From Wedeven, et al. (ref. 4), usii^ the symbols of this report, the dimensionless 
inlet distance at the fully flooded - starved boundary can be written as 


m-^ = 1 + 3 . 52 



-,2/3 


H 


c,F 


(5) 


Comparing equation (3) with equation (5) indicates close agreement with reference 4. 
Wedeven, however, predicts a slightly higher value of the fully flooded - starved bound- 
ary than is predicted from the present results. 


FILM -THICKNESS FORMULAS 

Having clearly established the limitii^ location of the inlet boundary for the fully 
flooded condition (eqs. (3) and (4)), an equation defining the dimensionless central film 
thickness for lubricant starvation conditions will be developed. The relation between the 
dimensionless central film thickness in the starved and fully flooded conditions can be 
expressed in general form as 


H 

H 


c^ 

c,F 



( 6 ) 


Table ni shows how the ratio of the dimensionless inlet distance parameter to the 
fully flooded - starved boundary (m - l)/(m* - 1) affects the ratio of central film thick- 


6 



ness in the starved and fully flooded conditions F* ^ least-square power curve 

fit to the 16 pairs of data points ’ ’ 



where i = 1,2, , . . , 16 


was used in obtainit^ the dimensionless central film thickness for a starved condition 


H 


c,S 



(7) 


By using a similar approach while making use of the data in table III, the dimensionless 
minimum film thickness for a starved condition can be written as 


H = H 

min,S “ min 




0.25 


( 8 ) 


Therefore, whenever m < m* , where m* is defined by either equation (3) or (4), a 
lubricant starvation condition exists. When this is true, the dimensionless central film 
thickness is expressed by equation (7) and the dimensionless minimum film thickness is 
expressed by equation (8) . If m2: m* , where m* is defined by either equation (3) 
or (4), a fully flooded condition exists. The dimensionless central and minimum film 
thicknesses for a fully flooded condition (H^ p and p) were developed in refer- 
ence 8 and can be written as ’ ’ 


H 


■P = 2.69 - 0.61 

C, r 

H Q ftQ ttO* 68|-(0 • 49ttt— 0 . 073/^ — 0.68k\ 

Hmin,F = 3.63 U G W (1 - e ) 


(9) 

( 10 ) 


where the dimensionless speed parameter 


U = 


E’R 

X 


( 11 ) 


the dimensionless load parameter 
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the dimensionless material parameter 

Piv,as 

the dimensionless ellipticity parameter 

k = i 
b 

Equations (13) and (14) can be expressed in a more usable form as 

G = E’a 


( 12 ) 


(13) 


(14) 


(15) 



In equation (16) the ellipticity parameter is expressed strictly in terms of the radii of 
curvature and thereby eliminates the common practice of evaluating the elliptical inte- 
grals of the first and second kinds. 

The ratio of the dimensionless inlet distance to the fully flooded - starved boundary 
as obtained from Wedeven, et al. (ref. 4), expressed as (m - l)/(m^ - 1), is also given 
in table in. By comparing these results with the results obtained from the present work, 
(m - l)/(m* - 1), it can be seen that for group 1 the agreement is excellent. However, 
the agreement in groups 2 and 3 is not as good. A possible explanation for this differ- 
ence can be that in the Wedeven, et al. (ref. 4) analysis, an approximate expression for 
the Hertzian deformation is used. They indicate that their equation (eq. (5)) is only valid 
for small m* or more specifically m* < 3. In group 2, m* = 3. 71; and in group 3, 
m* = 5.57. Since no such assumption is required in deriving equations (3) and (4), they 
would seem to be more appealing. 

Figure 2 shows the influence of inlet boundary parameter upon central film thickness 
for the Wedeven, et al. (ref. 4) results and for those obtained from the present work. 
From this figure it is observed that the Wedeven, et al. (ref. 4) results give slightly 
higher values of the central -film -thickness ratio of starved to fully flooded conditions 
than those obtained from the present results. 
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CONTOUR PLOTS OF RESULTS 


To explain more fully what happens in going from a fully flooded to a lubricant star- 
vation condition, figures 3 to 8 are presented. In figvires 3(a), (b), and (c) contour plots 
of dimensionless pressure (P = p/E’) are given for group 1 of table II and for dimension- 
less inlet distances of 4, 2, and 1.25, respectively. In these figures and the remainder 
of the contour plots to be presented the symbol + indicates the center of the Hertzian 
contact. It should be remembered that because of the way the coordinates are made di- 
mensionless, the actual contact ellipse becomes a Hertzian circle regardless of the el- 
lipticity parameter . The Hertzian contact circle is shown in each figure by asterisks . 

At the top of each figure the contour label and its corresponding value are given. The in- 
let region is to the left and the exit region is to the right. 

In figure 3 the contour values were kept constant. In figxires 3(a) and (b) the pres- 
sure spikes are shown with the B and C conto^lr labels, respectively; and in figure 3(c) 
no pressure spike occims. This implies that as the dimensionless inlet distance m de- 
creases (or as starvation occurs) the pressure spike is suppressed. In figure 3(a), with 
m = 4, a fvilly flooded condition exists. In figure 3(b), with m = 2, a starved condition 
exists. And in figure 3(c), with m = 1.25, the starvation is even more severe. That is, 
once starvation occurs the severity of the situation increases as m is decreased. Note 
in figure 3(c), the most severe starvation case, that the contour labeled H does not ex- 
tend as far to the left as it did for the fully flooded pressure results shown in figure 3(a) . 

In figures 4(a), (b), and (c) contour plots of dimensionless film thickness (H = h/R^^) 
are given for group 1 of table II and dimensionless inlet distances m of 4, 2, and 1.25, 
respectively. These film -thickness results correspond to the pressure results shown in 
figure 3. The central portion of the film -thickness contours has become more parallel 
as starvation has increased (figs. 4(a) to (c)), with the minimum film area decreasing. 
Note also that the film -thickness contour values for the starved condition (fig. 4(c)) are 
much lower than those for the fully flooded condition (fig. 4(a)) . 

In figures 5(a), (b), and (c) contour plots of dimensionless pressure (P = p/E’) are 
given for group 3 of table II and dimensionless inlet distances m of 6, 2. 5, and 1.5, 
respectively. This set of results should be compared with those presented for the lower 
speed parameter U in figure 3. The contour values in figures 5(a), (b), and (c) are the 
same in each case. In figure 5(a) a fully flooded condition exists, while in figure 5(c) a 
severely starved condition exists. The following observations can be made about fig- 
ure 5: 

(1) The distance from the center of the contact to the upstream location of the larg- 
est contour (labeled H) decreases as the severity of lubricant starvation increases. 

(2) Contours A, B, and C are absent in the severely starved condition (fig. 5(c)). 

Contour plots of the dimensionless film thickness (H = h/R^j.) are shown for group 3 

of table II and for dimensionless inlet distances m of 6, 2.5, and 1.5, respectively, in 
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figures 6(a), (b), and (c). These film -thickness results correspond to the pressure re- 
sults shown in figure 5. In figure 6(a) (the fully flooded condition) the minimum film 
area occiars to the sides of the conjunction in two areas that are midway between the cen- 
ter of the contact and the Hertzian circle . In figure 6(c) (the severely starved condition) 
the central portion has roughly parallel contours in the direction of motion with one min- 
imum area directly behind the axial center of the contact and near the Hertzian circle. 
Ranger, et al. (ref. 9) also found a similar contour plot in their 1975 publication. Note 
the similarity between the film -thickness contours of figure 4(c) (group 1 of table II) and 
that of figiire 6(c) (group 3 of table II) , From the labels of the contours in figure 6 we 
see that the film thickness for the starved condition is much lower than the film thick- 
ness in the fully flooded condition. 

The dimensionless pressure (P = p/E’) on the X-axis is shown for three values of 
dimensionless inlet distance m and for groups 1 and 3 of table II, respectively, in fig- 
iires 7(a) and (b). The value of Y is held constant near the axis of symmetry of the 
conjunction. In these figures as a conjunction becomes starved (as m is decreased) the 
pressvure spike diminishes. 

Figures 8(a) and (b) show the dimensionless film thickness (H = h/R^^) on the X-axis 
for three values of the dimensionless inlet distance m and for groups 1 and 3 of table H, 
respectively. The value of Y is held fixed close to the axis of symmetry of the contact. 
In these figures, particularly figure 8(b), the central region becomes flatter as starva- 
tion occurs. Also, in going from a fully flooded condition to a starved condition the film 
thickness decreases substantially. 


CONCLUDING REMARKS 

By using the theory and numerical procedure outlined and developed by the authors 
in earlier publications, the influence of lubricant starvation upon minimum film thickness 
in starved elliptical elastohydrodynamic conjunctions has been investigated. This study 
of lubricant starvation was performed simply by moving the inlet boimdary closer to the 
center of the conjunction. From the restilts it was found that the location of the dimen- 
sionless inlet boundary m* between the fully flooded and starved conditions could be ex- 
pressed simply as 


m* 



0.58 


H 


C,F 


or 
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where 

H p. fully flooded dimensionless central film thickness 
c, r 

Hmin F fully flooded dimensionless minimum film thickness 
m dimensionless inlet distance 

m* dimensionless inlet distance at the fully flooded - starved boundary 

Contour plots of the pressure and the film thickness in and around the contact have 
been presented for both the fully flooded and starved lubrication conditions. It is evident 
from the contour plots that the pressure spike becomes suppressed and the film thickness 
decreases substantially as the severity of starvation increases. 

The results presented in this report, when combined with the findings previously re- 
ported, enable the essential features of starved, elliptical elastohydrodynamic conjunc- 
tions to be ascertained. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 28, 1976, 

505-04. 
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TABLE I. - EFFECT OF STARVATION ON MINIMUM 


FILM THICKNESS FOR THREE GROUPS OF THE 
DIMENSIONLESS SPEED AND LOAD PARAMETERS 


Dimension- 


Group 


less inlet 
distance, 

1 

2 

3 

m 

Dimensionless load parameter, W 


0.3686x10'® 

0.7371x10'® 

0.7371X10*® 


Dimensionless speed parameter, U 


0.1683x10'^^ 

1.683x10'^^ 

5.050x10*^^ 


Minimum film thickness, 

’ min 

6 

6.317x10*® 

29.75x10'® 

61.32x10'® 

4 

29.27 

57.50 

3 

6.261 

27.84 

51.70 

2.5 


26.38 

46.89 

2 

5.997 

23.46 

39.91 

1.75 


21.02 

34.61 

1.5 

5,236 


27.90 

1.25 

3.945 




TABLE II. - EFFECT OF DIMENSIONLESS SPEED AND LOAD PARAMETERS ON 
DIMENSIONLESS INLET DISTANCE AT FULLY FLOODED - STARVED BOUNDARY 


Group 

Dimension- 
less speed 
parameter, 
U 

Dimension- 
less load 
parameter, 
W 

Dimension- 
less radius 
parameter , 

R^A. 

1 

0. 1683x10"^^ 

0.3686x10’® 

205.9 

2 

1.683 

.7371 

163.5 

3 

5.050 

.7371 

163.5 


Fully flooded 
central film 
thickness, 

«c,F 

Fully flooded 
minimum 
film thick- 
ness, 

^min, F 

Dimension- 
less inlet 
distance at 
fully flooded - 
starved 
boundary , 
m* 

7.480x10’® 

6.221x10’® 

2.62 

33.55 

29.20 

3.71 

70.67 

60.92 

5.57 




TABLE in. - EFFECT OF DIMENSIOKLESS INLET DISTANCE ON DIMENSIONLESS 


CENTRAL AND MINIMUM FILM THICKNESS RATIOS 


Group 

Dimension- 

Ratio of central 

Ratio of minimum 

Inlet boundary 

Inlet boundary 


less inlet 

film thicknesses 

film thicknesses 

parameter, 

parameter of 


distance, 

for starved and 

for starved and 

(m - l)/(m* - 1) 

Wedeven, et al. 


m 

flooded condi- 

flooded conditions, 


(1971), 



tions, 

®c,s/®c,F 

®tnin,s/®min, F 


(m - l)/(in^ - 1) 


2.62 

1 

1 

1 

0.9895 


2 

.9430 

.9640 

.6173 

.6108 


1.5 

.7697 

.8417 

.3086 

.3054 


1.25 

.5689 

.6341 

.1543 

.1527 

2 

3.71 

1 

1 

1 

0.8281 


3 

.9574 

.9534 

.7380 

.6111 


2.5 

.8870 

.9034 

.5535 

.4584 
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Ratio of central film thickness for 
starved and fully flooded conditions. 


Edges of computing area, where pressure is assumed ambient.-* 
Figure 1. - Computing area in and around the Hertzian contact circle. 



Figure 2. - Influenceof inlet boundary parameter upon central-film 
thickness ratio. 
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Dimensionless 

pressure, 

P - p/E' 

A 1.17x10'^ 

B 1.14 
C 1.10 
D 1.00 
E .90 
F .70 
G .40 
H .10 














Dimensionless 
film thickness, 
H-h/Rx 

A 4.(K10‘^ 

B 4.2 
C 4.5 
D 4.9 
E 5.4 
F 6.0 
G 6.8 
H 7.8 



(c)m»1.25. 


Figure 4. - Concluded. 


















Dimensionless 
film thickness. 
H * h/Rj( 

A 28.0x10"^ 

B 28.5 
C 29.0 
D 30.0 
E 32.0 
F 35.0 
G 39.0 
H 44.0 



(c) m “ 1.5. 


Rgure 6. - Concluded. 
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(a) Dimensionless parameters U and W heW constant as given in 
group 1 of tabie II. 

Figure 7. - Dimensionless pressure on X-a^s for three values of di- 
mensionless inlet distance. The value of Y Is held fixed near axial 
center of contact. 
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(b) Dimensionless parameters U and W held constant as given in 
group 3 of table II. 

Figure 8. - Concluded. 
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